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Work  on  this  new  program  was  effectively  initiated  in  April/May,  and  the  body  of  this  report 
will  mainly  be  concerned  with  new  initiatives  and  completion  of  earlier  ONR  funded  work. 

While  no  definite  condlusions  have  yet  been  reached,  the  basic  nature  of  the  problems  to  be 
studied  has  been  defined. 

A.  Goals 

The  objective  of  this  new  program  can  be  seated  very  simply;  it  is  to  incorporate  the  effects  of 
electronic  charge  redistribution  into  calculations  of  changes  in  lattice  potential  energies  for  arbitrary 
crystal  distortion  and  thus  to  predict  its  influence  on  the  uniform  polar  (antiphase  motion  of  positive 
and  negative  ions)  distortions  associated  with  ferroelectricty  in  oxide  based  ferroelectrics. 

Our  belief  that  this  is  central  to  the  understanding  of  ferroelectric  behavior  in  these,  the  over¬ 
whelming  majority  of  such  materials,  is  based  on  earlier  studies  during  the  previous  program. 
These  demonstrated  two  basic  facts:  (a)  Theoretical  Studies  on  halide  based  perovskites  revealed 
that  some  of  these;  e.g.,  NaCaF3,  can  show  ferrolectric  behavior  purely  due  to  ionic  “size”  effects 
(electronic  polarization  is  secondary  due  to  the  small  polarizability  of  the  F“  ion).  However,  this 
transformation  occurs  because  the  Na+  and  Ca2+  are  approximately  equal  in  size  and,  unfortu¬ 
nately,  this  condition  apparently  precludes  the  existence  of  NaCaF3  as  an  equilibrium  phase  under 
ambient  pressure,  (b)  Similar  theoretical  investigations  of  BaTiC>3  revealed  no  ferroelectric  insta¬ 
bility:  certainly,  the  widely  disparate  sizes  of  Ba2+  and  Ti4+  rules  out  behavior  similar  to  that  of 
NaCaFs.  Furthermore,  the  most  obvious  other  possibility,  that  Ti4+  is  so  small  that  it  “rattles” 
inside  its  octahedral  “cage”  of  six  oxygens  is  again  ruled  out  since  the  ionic  radii  are  such  that 
Ti4+  and  02-  come  into  contact  marginally  before  adjacent  pairs  of  O2-  ions  overlap.  Again,  one 
suspects  that  attempts  to  reverse  this  by  using  different  ions  could  well  render  the  perovskite  phase 
thermodynamically  unstable. 

Evidently,  there  is  something  missing  in  the  description  of  the  oxide  based  systems  and  the 
obvious  candidate  is  electron  redistribution. 

However,  it  has  been  known  for  ~  50  years  that  attempts  to  approximate  this  by  point  dipoles 
at  the  ion  centers  produces  completely  unphysical  results,  even  in  the  simplest  systems,  e.g.,  the 
heavier  sodium  halides  are  no  longer  stable  in  the  rocksalt  structure.  Furthermore,  one  also  has  to 
explain  why  BaTiC>3  is  strongly  ferroelectric  while  SrTiC>3  is  only  on  the  verge  of  being  polar. 

Very  recent  NRL  work  based  on  LAPW  total  energy  calculations  does  predict  a  ferroelectric 
ground  state  for  BaTi03  but  has  not  clearly  defined  the  origin  of  this  behavior,  apart  from  indi¬ 
cating  that  charge  transfer  from  the  O2-  ions  into  the  unfilled  d  states  of  Ti4+  plays  a  significant 


rolth-  Unfortunately,  such  studies  are  extremely  time  consuming,  and  offer  no  hope  of  providing 
the  type  of  potential  function  that  we  have  obtained  from  Gordon-Kim  rigid  ion  calculations  for 
halide  based  systems.  Clearly  incorporation  of  electron  redistribution  effects  effectively  increases 
the  number  of  degrees  of  freedom  per  ion,  but  this  can  be  handled  if  this  redistribution  can  be 
characterized  by  a  relatively  small  number  of  parameters.  An  example,  is  the  point  dipole  approx¬ 
imation,  alluded  to  earlier.  Evidently,  this  is  a  very  bad  approximation;  but  it  does  completely 
characterize  the  charge  redistribution  using  only  three  additional  degrees  of  freedom  per  ion.  It 
is  clearly  implausible  that  any  description  can  have  fewer  parameters,  but  it  does  not  necessarily 
follow  that  any  triplet  of  parameters  will  be  equally  bad. 

Currently,  a  sequence  of  LAPW  calculations  has  been  initiated  at  NRL  designed  to  examine 
in  detail  the  nature  of  the  charge  distortions  in  BaTiC>3  associated  with  polar  motion.  The  first 
objective  is  to  determine  the  best  representation  of  the  charge  distribution  in  terms  still  of  purely 
spherical  ions.  Then  the  next  stage  will  be  to  isolate  the  regions  in  which  the  fit  is  worst  and  then 
to  further  reduce  the  discrepancy  by  introducing  new  spherical  distributions  away  from  the  ion 
centers. 

B.  Results 
Initiation  of  New  Work 

The  optimized  spherical  distributions  alluded  to  in  the  previous  section  were  derived  at  NRL 
using  the  “Potential  Induced  Breathing”  (PIB)  model.  As  a  first  stage  in  the  present  study, 
Professor  J.  W.  Flocken  has  been  applying  this  to  LaA103  in  order  to  see  if  we  can  correctly 
reproduce  the  cell  multiplying  transition  in  that  system.  Results  thus  far  indicate  considerable 
sensitivity  to  the  manner  in  which  the  02-  ion  is  stabilized  in  PIB. 

Dr.  H.  M.  Lu  has  been  studying  the  dynamics  of  ionic  molecular  systems  when  one  simulates 
polarization  effects  by  permitting  rigid  displacement  of  the  valence  electrons  with  respect  to  the 
core  but  uses  a  priori  potentials  for  the  various  interactions.  This  work  is  ongoing  and  could  well 
be  generalized  to  introduce  other  spherical  charge  distributions  not  centered  on  the  ions; 

Professor  R.  A.  Guenther  has  begun  extending  to  two  dimensions  the  theory  of  the  dielectric 
response  of  double- well  systems,  which  was  central  to  our  understanding  of  dielectric  loss  in  SbSI 
under  the  previous  Award. 

Professor  J.  W.  Flocken  and  Dr.  Z  Mo  have  been  completing  their  studies  on  mixed  halide 
perovskites;  specifically,  the  (Na,  K)  CaF3  system.  Since  (non-ferroelectric)  KC&F3  can  be  grown, 


we  wish  to  discover  what  level  of  Na  substitution  is  needed  to  initiate  the  ferroelectric  behavior 
predicted  theoretically  for  NaCaFs  since,  if  this  is  low  enough,  it  should  be  possible  to  produce 
such  systems  by  doping  pure  KCaF3. 

The  three  graduate  students  who  have  begun  work  since  the  initiation  of  the  new  program  are 
principally  acquiring  familiarity  with  our  existing  codes  and  molecular  dynamics  capabilities.  They 
are  also  learning  the  basic  theory  underlying  these  tools.  This  provides  an  opportunity  for  them  to 
begin  exploratory  studies  on  the  general  problem  of  obtaining  accurate  representation  of  electron 
charge  distribution  with  a  minimal  number  of  parameters. 

C.  Future  Work 

During  the  coming  year,  work  on  the  topics  just  described  will  be  continued  with  particular 
emphasis  uniting  the  NRL  studies  on  LAPW  charge  densities  with  our  work  on  representing  these 
densities  and  variations  with  minimal  numbers  of  parameters  determined  by  self-ccnsistency.  It  is 
possible  that  other  topics  may  emerge  as  priorities  thus  dictating  new  areas  of  emphasis;  however 
in  view  of  the  limited  duration  of  the  present  program  to  date,  it  would  be  unwise  at  this  stage  to 
instance  specific  possibilities. 
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